Gravothermal Collapse of Self-Interacting Dark Matter Halos 
and the Origin of Massive Black Holes 
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Black hole formation is an inevitable consequence of relativistic core collapse following the gravother- 
mal catastrophe in self-interacting dark matter (SIDM) halos. Very massive SIDM halos form super- 
massive black holes (SMBHs) > 10® Mq directly. Smaller halos believed to form by redshift 2 = 5 
produce seed black holes of lO'^ — 10^ Mq which can merge and/or accrete to reach the observational 
SMBH range. This scenario for SMBH formation requires no baryons, no prior star formation and 
no other black hole seed mechanism. 

Pacs Numbers: 95.35.-fd, 97.60.Lf, 98.62. Ai, 98.62. Gq 

Black holes in the centers of galaxies are inferred to 
have masses of 10^ — 10^ Mq, and to date dozens of 
candidates have been discovered Black holes in this 
mass range are the likely power sources in quasars and 
active galactic nuclei Recent observational data 

have identified a strong correlation between the inferred 
black hole mass and the stellar velocity dispersion in the 
host galaxy bulge a correlation has also been found 
between the black hole mass and the mass of the bulge . 
Both correlations suggest that formation and evolution of 
the central black hole and the bulge of the host galaxy 
may be closely related. 

The origin of these supermassive black holes (SMBHs) 
is uncertain. Proposed scenarios include the collapse of 
a supermassive star 0, possibly built up by collisions 
and mergers of ordinary stars Another scenario in- 
volves the collapse of the core of a dense relativistic clus- 
ter (e.g. neutron stars or stellar mass black holes) follow- 
ing gravothermal evolution. The cluster core evaporates 
mass to an extended halo on a gravitational scattering 
relaxation timescale, while the core density and velocity 
dispersion grow (the "gravothermal catastrophe" ; P , p^ ). 
Zel'dovich and Podurets originally conjectured pd[ ] and 
Shapiro and Teukolsky subsequently demonstrated 
that a nearly coUisionless gas in virial equilibrium like a 
star cluster experiences a radial instability to collapse on 
a dynamical timescale when its core becomes sufficiently 
relativistic. As the instability sets in, the core and its im- 
mediate surroundings undergo catastrophic collapse to a 
black hole and the ambient halo settles into a new dy- 
namical equilibrium state around the central hole. 

Here we report that the formation of a black hole in the 
center of a galactic halo is a natural and inevitable conse- 
quence of the gravothermal evolution of self-interacting 
dark matter (SIDM) halos. The possibility that dark 
matter particles are "self-interacting" has been revived 
recently Studies of SIDM via N-body simulations 

p^ , p^ and via a gravothermal approach have con- 
firmed that SIDM halos are more consistent with observa- 
tions: they exhibit a flat density core, rather than a cuspy 



one that arises for cold dark matter [ |l7yi§| ] . An isolated 
SIDM halo evolves gravothermally, since the thermal re- 
laxation timescale due to collisions is shorter than the 
typical ages of cosmological halos. Unlike a star clus- 
ter, a SIDM halo retains an appreciable core mass as it 
evolves towards the relativistic instability Q . 



I. GRAVOTHERMAL EVOLUTION AND 
COLLAPSE OF A SIDM HALO 

The structure of a relaxed SIDM halo can be defined 
by the central mass density pc, the (one-dimensional) ve- 
locity dispersion v^, and the total halo mass Mtot- For 
a given cosmology, the background density, p, is a func- 
tion only of redshift. We assume that the central density 
contrast x = pc{t = 0)/p at halo virialization is indepen- 
dent of redshift and halo mass, which is consistent with 
the simulations of Dave et al. |jl^ and recent work on 
truncated isothermal spheres by Shapiro, Iliev and Raga 
]l9[ . The independent parameters of a SIDM halo are 
then Mtot, x and zq, the redshift of virialization. Define 
mass and radius scales, Mq and Rq respectively, where 
Mo = 4:TtRqPc and = GMq/Rq, and also the ratio 
g = Mtot/Afo, which is a unique function of x. The cen- 
tral velocity satisfies 

v^t = 0) = (4^)i/=^G(Mtot/5(2:))'/=^(xp(zo))'/' . (1) 

Evolutionary Stages. The properties in the halo de- 
pend critically on the ratio of the coUisional mean-free- 
path, A, to the gravitational scale height, H: 



A = (pa)- 
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where p is the local mass density, v is the local one dimen- 
sional velocity dispersion, and a is the SIDM cross section 
per unit mass. In the long mean-free-path (Imfp) region 
A 3> -ff, particles perform several orbits between colli- 
sions; in the short mean-free-path (smfp) region A <C -ff , 
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particle motion is severely restrained, and heat transfer 
proceeds through diffusion, as in a fluid. Simulations 
show that a halo core must form in the Imfp regime ; 
otherwise the particles generate strong shocks, whereby 
entropy loss to the halo destroys the central core struc- 
ture @. 

The requirement that at formation the core satisfies 
{X/H)c > 1 yields an upper limit for the halo mass: 



Mtot < 47r(/9(zo)a;) ^gix)a' 



(3) 



Black Hole Formation Gravothermal evolution of a 
SIDM halo is shown in Figure While the extended 
halo {p oc r^^-^^) remains in the Imfp limit, the core be- 
comes increasingly dense. The core eventually bifurcates 
into a smfp inner part which is fluid- like (A <^ H), and 
a Imfp outer part which is nearly static. The transition 
region corresponds to A/iJ ~ 1. It is the inner core and 
its immediate surroundings which collapse to a black hole 
following the onset of the relativistic instability, leaving 
the outer core and extended halo in dynamical equilib- 
rium about the central hole fl^ . 

Evaporation of mass due to collisions from a Imfp core 
into the halo is rapid with d\ogMc/d\og{v^) « —4.27 
. But once the core becomes very dense, with central 
X/H > 100, evaporation is limited to a surface effect, 
which is much less efficient. We find that if the core 
of a SIDM halo forms with X/H > 1, its mass decreases 
by about one order of magnitude prior to reaching the 
smfp limit (while the central velocity hardly changes). 
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FIG. 1. Snapshots of the (a) density and (b) velocity dis- 
persion profiles of a SIDM halo with a = 0.6747r7?o/M'o 
(for a — 2.26) at selected times during gravother- 
mal evolution [16]. The thick line denotes the pro- 
file at t = 0. Subsequent profiles correspond to 
t/tr{0) =264.0, 287, 289.5, 289.90, 290.12, 290.30, 290.49, and 
290.63. 



after which the core mass decreases according to 
dlogMc/cilog(w^) ~ —0.85. This estimate gives the in- 
ner core mass Mcou at the onset of relativistic instability 
(vc « c/3): 
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Core Lifetime. The lifetime of the core until collapse is 
Jl6t tcoii ~ 290ir(0), where tr(0) is the collision relaxation 
timescale in the core at formation. 
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(pc(i = 0)^;c(t = 0)a)" 



(5) 



and a ~ 2.26 for purely elastic, hard-sphere collisions and 
a Maxwell-Boltzmann velocity distribution. Requiring 
that the core of a halo formed at redshift zq collapses by 
redshift zi yields a lower limit on the mass of the SIDM 
halo: 



m,„,>(^V(4.g3)-/.M£M:!{^ 

a J ^ ' [t{z,)-t{zo)f 



(6) 



II. A DIRECT SCENARIO FOR SMBH 
FORMATION 

Since the core virializes with nonrelativistic veloci- 
ties, Vc < 10^ — lO'^ km s~^ , the core mass at col- 
lapse will be several orders of magnitude smaller than 
Mc(0) w 10~2Mtot according to Eq. (||). Consequently, 
gravothermal evolution directly culminating in a SMBH 
with a mass above 10^ Mq requires that the total mass 
of the progenitor SIDM halo exceed 10^^ Equa- 
tions (^) and (^) determine the range of halo masses 
formed at redshift zq that undergo core collapse by red- 
shift zi. These masses increase for smaller values of the 
ratio X, i.e., for smaller central density contrasts at halo 
formation. 

We examine this scenario by setting x to the relatively 
small value of a; = 1.8 x lO'', which leads to g « 206. This 
is the value derived for a truncated isothermal sphere of 
dark matter |l9) for an Einstein de-Sitter universe, and 
also applies for all but very small redshifts for a flat uni- 
verse with a finite cosmological constant |^l| . In Figure ^ 
we show the range of halo masses which experience core 
collapse (upper frame) as a function of formation redshift 
Zq. Since the lower limit depends also on the desired red- 
shift at collapse, zi, we show three variations: zi = 6,3 
and 0. The lower frame shows the corresponding core 
(« black hole) masses at collapse, using Eqs. (1) and 
(^). The cosmological model is flra = 0.3, f^A — 0.7 and 
h — 0.65, and we set a = 5 cm^ gm^^ p^]. 
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FIG. 2. Masses of SIDM halos that can undergo core col- 
lapse by redshifts zi — 6, 3 and (upper frame), and the 
corresponding core mass at collapse (lower frame), as func- 
tions of formation redshift zq- Dashed curve corresponds to 
the upper limit (Eq. [^) and solid curves to the lower limits 
(Eq. Iq]) as functions of zq and zi. In this model the density 
contrast x = pc(t — Q)/p— 1.8 x 10*. 

We find that SMBHs of 10*^ - 10^ Mq can form in mas- 
sive halos vifhich experience core collapse at intermediate 
and low redshifts. At redshift zi = 6, our model limits 
the initial black hole mass to M w 5 x IO^A/q. Subse- 
quent accretion onto the black hole must be invoked to 
reach masses of 10^ Mq |22|, especially to reconcile with 
recent observations of high red shift quasars |Q . Very 
young halos are practically excluded as candidates for 
reaching core collapse by redshift z — 0: halos formed 
later than z ~ 1.5 will not have reached core collapse, 
regardless of mass, and only cores of extremely massive 
halos (Aftot > 10^^ Mq) could have collapsed if formed 
after z — ?>. 




FIG. 3. Same as Figure ^ except that 
a; = pc(t = 0)/p= 10^ 

The relations Afmin oc a;"^/^ and Af,„ax oi x"^ show 
why core collapse shifts to relatively low mass halos for 
X = 10^. For this x any halo which virializes at redshift 
z > 5 will have undergone core collapse. However, young 
halos formed at low redshifts will still be safe from core 
collapse, as we found also for x = 1.8 x 10'*. In this 
scenario, core collapse even by redshift z 6 is quite 
possible, for halos with masses in the range 10^ — 10*° Mq 
virialized at redshifts z = 7 — 10. The masses of the 
collapsed cores are small, and for halo masses considered 



coll 



< 



realistic in the Press-Schechter formalism [ 
10^ — 10^ Mq. In this case gravothermal core collapse 
does not lead directly to the formation of SMBHs found 
in the centers of galaxies, but the resulting intermediate 
mass black holes can serve as seeds for SMBH build-up 
through mergers 22|, or accretion. 



III. A BOTTOM-UP SCENARIO FOR SMBH 
FORMATION 

Large values of Afcoii require very massive halos in the 
direct scenario. Virialization of massive halos, especially 
at high redshift, is not favored in the currently accepted 
picture of structure formation ||2j,|2^ . The smaller mass 
halos formed preferentially in current theory can undergo 
core collapse in a Hubble time provided the fiducial value 
of a; is larger (see Eqs. PH). In Figure ^jwe show similar 
results to those of Fig. y, except that now x = 10® (and 
g « 902). This value appears to be in better agreement 
with the numerical simulations of ref. at least for 
low redshifts. 



IV. DISCUSSION 

Gravothermal core collapse in a SIDM halo triggers the 
formation of a SMBH at its center. The black hole forms 
when the inner core becomes relativistic and dynamically 
unstable. The mass of the collapsing core will be 10^* — 
10^® times the total mass of the halo. Forming SMBHs 
by core collapse in SIDM halos requires no baryons, no 
prior epoch of star formation and no other black hole 
seed mechanism. 

For massive (> 10^^ Mq) halos which virialize at mod- 
est overdensities, SMBHs with masses of 10® — 10^ Mq 
form directly through gravothermal collapse. An alterna- 
tive scenario where the smaller SIDM halos preferred in 
current cosmological models reach core collapse in a Hub- 
ble time is also possible and arises if the typical overden- 
sity of the halos at virialization is large. In this case any 
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halo which forms prior to redshift z = 5 produces a black 
hole by z = 0, and black hole formation at high redshift 
is also possible. However, these halos give rise only to 
low and intermediate mass black holes, < 10^ — 10^ Mq. 
These might be seeds of SMBHs through multiple halo 
mergers |2^. Even if only a few percent of the halos 
achieve core collapse by z — 5, this initial black hole 
population is sufhcient to generate the observed SMBH 
spectrum eventually ||2^ . Coalescence of black holes in 
halo mergers is a likely source of gravitational waves, po- 
tentially detectable with the Laser Interferometer Gravi- 
tational Wave Observatory (LIGO) as well as the planned 
Laser Inferometer Space Antenna (LISA). If most of the 
SMBH population arises through multiple mergers of 
lower mass seed black holes, the rate may be as large 
as several per year p6| . 

The newly formed black hole may grow through accre- 
tion of SIDM from the halo. The main reservoir for accre- 
tion is the outer core, which is characterized by A/iJ > 1, 
and includes a mass ^ 10~'^Mtot. If the ambient SIDM 
halo eventually dominates the spherical bulge of galaxies, 
then the accretion of the entire outer core would provide 
a natural means of producing the ratio MenZ-^^buige ob- 
served in nearby galaxies. It would also explain the origin 
of the most massive (> 10^ Mq) central black holes. The 
possibility that SIDM accretion onto seed black holes is 



the origin of SMBHs was also proposed in |27|, but there 
it is assumed that the seeds arise from supernovae ex- 
plosions in massive stars rather than the gravothermal 
scenario discussed here. Their quantitative results de- 
pend crucially on an assumed singular power-law density 
profile of the SIDM, which has yet to be verified self- 
consistently. 

Young halos, especially low mass ones, could not have 
reached core collapse by z = 0. This allows SIDM ha- 
los to explain the flat density cores observed in some 
dwarf and low-surface- brightness galaxies, where the in- 
ferred values of central density and velocity dispersion 
are pc ~ 0.02 M© pc^^ and Vc ^ 10'' cm s~^ , yield- 
ing a gravothermal core lifetime that greatly exceeds the 
Hubble time. 

To identify the main route for SMBH formation, it 
would be useful to determine the SMBH population as 
a function of redshift. More data about the presence or 
absence of central black holes in dwarf galaxies would 
also be useful. 
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